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Abstract 

We consider the semiclassical limit for the Heisenberg-von Neumann equation with a 
potential which consists of the sum of a repulsive Coulomb potential, plus a Lipschitz po- 
tential whose gradient belongs to BV; this assumption on the potential guarantees the well 
posedness of the Liouville equation in the space of bounded integrable solutions. We find 
sufficient conditions on the initial data to ensure that the quantum dynamics converges to 
the classical one. More precisely, we consider the Husimi functions of the solution of the 
Heisenberg-von Neumann equation, and under suitable assumptions on the initial data we 
prove that they converge, as e — > 0, to the unique bounded solution of the Liouville equation 
(locally uniformly in time). 

1 Introduction 

The aim of this paper is to study the semiclassical limit for the Heisenberg-von Neumann (quan- 
tum Liouville) equation: 

' iedtpl = [H e ,fi,], 

(1.1) 

, P0 = P0,e, 

{po,e}e>o being a family of uniformly bounded (with respect to e), positive, trace class operators, 
and with H £ = -^A + U. 

When p~o )£ is the orthogonal projector onto ^o,e £ L 2 (M. n ), (|1.2j) is equivalent (up to a global 
phase) to the Schrodinger equation 

(1.2) 

We recall that the Wigner transform W e ip of a function ip S L 2 (M. n ) is defined as 

W e ^{x,p) := — -!— / ij){x + £ -y)^{x - £ -y)e- ivv dy, 
(2vr) n J R n 2 2 
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and the one of a density matrix p is defined as 



W £ p(x,p) : 



(27T)" J w 



•71 



p(x + -y,x- -y)e %vv dy 



(1.3) 



where p(x, x') denotes the integral kernel associated to the operator p. 

The weak limit of the Wigner function of the solution of (|1.2p or (jl.ip has been studied in 
many articles (e.g. |15| ITS"! 114] . and more recently in strong topology in j6j |7]). More precisely, 
it is well-known that the limit dynamics of the Schrodinger equation is related to the Liouville 
equation 



of distribution to some (nonnegative) measure po, then W £k p £ t k — > (3>t)#/io (the convergence 
is again in the sense of distribution), where <£>£ is the (unique) flow map associated to the 
Hamiltonian system 



so that pt '■= i s the unique solution to fj 1 . 4 [) (here and in the sequel, # denotes the 

push-forward, so that p t (A) = p (^ 1 (A)) for all A C R 2n Borel). 

(B) If U is of class C 1 and there exists a sequence Eu — > such that the curve t i— >■ W £k p £h 
converges in the sense of distribution to some curve of (nonnegative) measure t i— > pt, then pt 
solves ()1.4jl . 

In the present paper we want to use some recent results proved in |31 H] to improve the 
literature in two directions: 

(i) By lowering the regularity assumptions of (A) on the potential in order get convergence 
results for a more general class of potentials, as described below. 

(ii) Get rid of the "after an extraction of a subsequence" argument, due to compactness, used 
in most of the available proofs where one is unable to uniquely identify the limit. More 
precisely, in (B) above one needs to take a subsequence along which the whole curve 
1 1 — y W £k p £ t k converges for all t in order to obtain a solution to f j 1 . 4 [) . Moreover, the limiting 
solution may depend on the particular subsequence. In our case we will be able to show 
that, for a class potential much larger than C 2 , once one assumes that the Wigner functions 
at time t = have a limit, then the limit at any other time will converge to a "uniquely 
identified" solution of f j 1 . 4 1) . 

The price to pay for the lack of regularity of the potential will be to have some size condition on 
the initial datum which forbids the possibility of considering pure states. Even more, the Wigner 
function of the initial datum cannot concentrate at a point, a possibility which might actually 
enter in conflict with the fact that the underlying flow is not uniquely defined everywhere. Let us 
mention however that, with extra assumptions on the potential (but still allowing the possibility 
of not having uniqueness of a classical flow), it is possible to consider concentrating initial Wigner 





(1.5) 
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functions, giving rise to atomic measures whose evolution follows the "multicharacteristics" of 
the flow (see [7]). 

As described below, we will nevertheless show that, for general bounded and globally Lipschitz 
potential associated to locally BV vector fields (in addition to some Coulomb part), the Wigner 
measure of the solution at any time is the push-forward of the initial one by the Ambrosio- 
DiPerna-Lions flow [9j[2]. 

Our method will use extensively the Husimi transforms ip i— > W £ ip and p i— > W e p, which 
we recall are defined in terms of convolution of the Wigner transform with the 2n-dimensional 
Gaussian kernel with variance e/2: 

W £ ^:={W e ^)*Gf n \ W £ p:=(W £ p)*G? n \ G? n \x,p) := & ^ + '*' - = G^{x)Gt\p). 

{■Key 1 

(1.6) 

Of course, the asymptotic behaviour of the Wigner and Husimi transform is the same in the limit 
e — > 0. However, one of the main advantages of the Husimi transform is that it is nonnegative 
(see Appendix). 

Let us observe that, thanks to (jA.8[) . the L°°-norm of W £ ip can be estimated using the 
Cauchy-Schwarz inequality: 



However, this estimate blows up as e — > 0. On the other hand we will prove that, by averaging 
the initial condition with respect to translations, we can get a uniform estimate as e — > (Section 
13. 2p . This gives us, for instance, an important family of initial data to which our result and the 
ones in [T] apply (see also the other examples in Section [3| . 

2 The main results 
2.1 Setting 

We are concerned with the derivation of classical mechanics from quantum mechanics, corre- 
sponding to the study of the asymptotic behaviour of solutions p\ to the Heisenberg-von Neu- 
mann equation 

iedtpf = [H £ ,p £ t ] 

(2.1) 

PO = P0,e, 

2 

as e — > 0, where H £ = — 4j-A + U, and U : R n — > R is of the form + U s on R n , where U s is a 
repulsive Coulomb potential 

U s (x)= V , ZiZj , , M < n/3, x = (xi,...x M ,x) G (M 3 ) M x R n ~ 3M , Z { > 0, 

' -J T ■ rp ■ 

l<i<j<M 1 Jl 
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Ub is globally bounded, locally Lipschitz, VUb <E BVi oc (R n ; M n ), and 

\VU b (x) 



ess sup 

xei" 1 + \ x \ 



< +00. 



The formal solution of (j 1 . 2 1) is pf, where 



and its kernel is p|. Moreover, as shown for instance in [15], W e pf solves in the sense of distri- 
butions the equation 

d t W £P e t + P ■ V x W £ p £ t = g £ (U,p £ t ), (2.2) 



where £ £ {U,p) is given by 
£ e (U,p)(x,p) :-- 



(2vr)" J Rn 



U(x + ^y)-U{x-^y) 



p{ x+ £ -y^ x - £ -y)e-^dy. (2.3) 



Adding and subtracting VU(x) ■ y in the term in square brackets and using ye tp ' y = iV p e tp ' y , 
an integration by parts gives (o £ (JJ, p) = VC/(x) • VpWsp + ^(U, p), where ^(U, p) is given by 



4(U,p)(x,p) :-- 



U{x + ^y)-U{x-^y) 1 e e . 

, VU{x) -y p{x + -y,x- -y)e py dy. 

(2.4) 

Let b : M. 2n — > M? n be the autonomous divergence-free vector field b(x,p) := (p, — VC/(x)). Then, 
by the discussion above, W e pf solves the Liouville equation associated to b with an error term: 



dtWefi + b- VW e p £ t =4(U, pt). 



(2.5) 



On the other hand, thanks to (|2.2p . it is not difficult to prove that W e p\ solves in the sense of 
distributions the equation 



where 



d t W £ p e t + p ■ V x W e p e t = g e (U, pt) * G^ - yfeV x ■ \W e p\ * Gf n \ 



(2.6) 



(2.7) 



Since W e pf and W e p e have the same limit points as e — > 0, the heuristic idea is that in the limit 
e — > all error terms should disappear, and we should be left with the Liouville equation (which 
describes the classical dynamics) 



dtUt + b ■ Vujt = on 



o'2n 
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2.2 Preliminary results on the Liouville equations 

Under the above assumptions on U one cannot hope for a general uniqueness result in the space 
of measures for the Liouville equation, as this would be equivalent to uniqueness for the ODE 
with vector field b (see for instance |3|). On the other hand, as shown in |T| Theorem 6.1], the 
equation 

dtOJt + b ■ Voj t = 

(2.8) 

lo = u) G L l (R 2n ) n L°°(R 2n ) and nonnegative, 

has existence and uniqueness in the space L^f([0, T\ \ L 1 (M 2n ) nL°°(IR 2 ™)). This means that there 
exist a unique W : [0, T] -> L 1 (M 2n )nL°°(]R 2n ), nonnegative and such that esssup 4g[0 T] ||^||li(r2„ ) + 
\\Wt ||L°°(R 2n ) < +oo, that solves (|2.8p in the sense of distributions on [0,T] x R 2n . 

One may wonder whether, in this general setting, solutions to the transport equation can still 
be described using the theory of characteristics. Even if in this case one cannot solve uniquely 
the ODE, one can still prove that there exists a unique flow map in the "Ambrosio-DiPerna-Lions 
sense". Let us recall the definition of Regular Lagrangian Flow (in short RLF) in the sense of 
Ambrosio-DiPerna-Lions: 

We say that a (continuous) family of maps &t '■ ^ 2n K 2n , t > 0, is a RLF associated to 
(L5]) if: 

- $o is the identity map. 

- For J2? 2n -a.e. (x,p), 1 1— > $t(x,p) is an absolutely continuous curve solving (11.51) . 

- For every T > there exists a constant C T such that ($t) # Jz? 2n < C T ^£ 2n for all t G [0, T], 
where J>f 2n denotes the Lebesgue measure on M. 2n . 

Observe that, since VU is not Lipschitz, a priori the ODE (j 1 . 5 [) could have more than one 
solution for some initial condition. However, the approach via RLFs allows to get rid of this 
problem by looking at solutions to (JT3J) as a whole, and under suitable assumptions on U the 
RLF associated to (j 1 . 5 [) exists, and it is unique in the following sense: assume that 4> x and <J> 2 
are two RLFs. Then, for Jz? 2n -a.e. (x,p), §\(x,p) = & 2 (x,p) for all t G [0,+oo). In particular, 
as shown in [I] Section 6], the unique solution to (|2.8p is given by 

c^ 2 ™ = (^) # (^ 2 "). (2.9) 

Hence, the idea is that, if we can ensure that any limit point of the Husimi transforms W £ pf 
give rives to a curve of measure belonging to L^°([0, T]; L 1 (]R 2n ) n L°°(IR 2n )), by the aforemen- 
tioned result we would deduce that the limit is unique (once the limit initial datum is fixed), 
and moreover it is transported by the unique RLF. In order to get such a result we need to make 
some assumptions on the initial data. 
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2.3 Assumptions on the initial data and main theorem 

Let {/Oo,e}ee(o,i) be a family of initial data which satisfy 

Po,e = Po, £ , Po,e>0 and tr(p , £ ) = 1 Vee (0,1). 

Let 

jeN 

be the spectral decomposition of po E , and denote by po,e its integral kernel. 
We assume: 

sup Y fifWH^W 2 <+oo, (2.10) 

^Po, e <CId, (2.11) 

lim sup / po e(^) a?) da; = 0, (2-12) 

fl->+oo £e (o,i) Jm.\b^ 

and 

where is the ball of radius R in W 1 and J 7 is the Fourier transform on M 2n , see (|A.5[) . 

Conditions (|2.12p and fj2 . 13[) are equivalent to asking that the family of probability measure 
{WePo,e}se(o,l) is tight (see Appendix). By Prokhorov's Theorem, this is equivalent to the com- 
pactness of {W e /9o,e} £ 6(o,i) with respect to the weak topology of probability measures (i.e., in 
the duality with Cb(M. 2n ), the space of bounded continuous functions). Hence, up to extract- 
ing a subsequence, assumptions (|2.12p together with (|2.13p is equivalent to the existence of a 
probability density co such that 

w - limW £ p 0i£ ^ 2n = u^ 2n £ ^(K 2n ), (2.14) 

where ^(]R 2n ) denotes the space of probability measure on M n . In order to avoid a tedious 
notation which would result by working with a subsequence e^, we will assume that f|2. 14[) holds 
along the whole sequence e — > 0, keeping in mind that all the arguments could be repeated with 
an arbitrary subsequence. 

Let us observe that condition (|2.10p is slightly weaker than sup eg ( 0i i) tr(H 2 p~Q^ £ ) < +oo, as 
in order to give a sense to the latter we need the operator H 2 p\ to make sense (at least on a 
core). Concerning assumption (|2.14p . let us observe that the hypothesis tr(po,e) = 1 implies that 
W e p , E G ^(R 2n ) (see Appendix). 

To express in a better and cleaner way the fact that the convergence is uniform in time, we 
denote by d^> any bounded distance inducing the weak topology in ^(R 2n ). Recall also that 
&t denotes the unique RLF associated to b(x,p) = (p, — V£7(x)), so that ($t)# (w,£? 2n ) is the 
unique nonnegative solution of fl23J) in L2>([0, T]; L 1 (IR 2n ) n L°°(lR 2n )). 



Theorem 2.1. Let U be as in Section HO , Under the assumptions &2.10\) . h2.11\) and {2-1$ 

limsupd^{W £ p £ t ^ 2n ,(^t)#(^ 2n )) = 0. (2.15) 
[o,r| 

Moreover, if we define WtJi? 2n = (<3?t)# (w.£? 2ra ) , for every smooth function ip G C^°(M 2n ) the 
map 1 1 — y J K 2„ f^i dx dp is continuously differentiate, and 



ipWt dxdp = b ■ VipWt dx dp. 



d 
dt 

The rest of the paper will be concerned with the proof of Theorem 12.11 However, before 
proceeding with the proof, we first provide some example and sufficient conditions for our result 
to apply. 



3 Examples 

We will give three types of examples of density matrices satisfying the assumptions of the pre- 
ceding section, so that Theorem 12. II applies. 



3.1 Average of an orthonormal basis 

For simplicity, we set up our first example in the one-dimensional case. In particular, there is 
no Coulomb interaction (that is, U = £/&), since by assumption Coulomb interactions are three- 
dimensional. We leave to the interested reader the extension to arbitrary dimension (the only 
difference in the case U s ^ appears when checking assumption (|2.10p ). 

Let us consider the orthonormal basis of L 2 (R) given by the (semiclassical) Hermite functions 

= fl ^ - Hj 3 € N, 

^ v / 2i7!(7re)V4 ] \^J 

where J9Vs are the Hermite polynomials, i.e. 

- l-V* ■ 

The following holds: 

Proposition 3.1. Let {/^ £ }yeN be a sequence of positive numbers, and define the density matrix 
p e given by 

Assume that 

• o< M f <ce, E i6 N^ 6) = 1 ; 
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Then (l2~T0l) . (I2TTT) . and (l2~T4l) /io/d. 

Proof. The first assumption is equivalent to (|2,lip and the trace-one condition. 
Concerning (|2,10p . using the well-know fact that 



by a simple calculation it follows that 

u e2 ^ /( £ ) , rr /( £ ) 

= -\ (VM^Vrfli - (2j + + v/(i + i)(i + 2)^ 2 ) + u b tf 

Hence 



I (Vj(j - 1) + (2j + 1) + v / (i + i)(i + 2)j + 



< C(l + e z j 

and (|2,10p follows from the first two assumptions. 

Finally, the third assumption implies (|2.14p by noticing that 

w - lim W E i/)f ) = 5(x 2 +p 2 -a) V a > 
(see, for instance, |15| Exemple III. 6]). □ 
3.2 Toplitz case 

Let cj) E H 2 (R n ; C) with / Rn \4>(x)\ 2 dx = 1. Given e, q > 0, for any </ G R n let be defined 
by 

Then, using Plancherel theorem, one can easily check that the identity 

\r w , q W w , g \dwdq = (27r)"Id (3.1) 

holds, where \ip){^\ is the Dirac notation for the orthogonal projection onto a normalized vector 
tp € L 2 (]R n ). Thanks to (|3.ip and the fact that orthogonal projectors are nonnegative operators, 



e e . 
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we immediately obtain the following important estimate: for every nonnegative bounded function 
Xe ■■ K 2n -> M, it holds 

Xs{w,q)W w>q )(i> £ w>q \ dwdq < || X£ ||oo(27r)"Id. (3.2) 
Set now 

Po,e-= Xe(w,q)\ip e w ){ilj e w \dwdq, e G (0, 1), 

where {Xe}ee(o,i) i s a family of nonnegative bounded functions such that J R2 „ Xe( w , q) dw dq = 1, 
and let S be the singular set of U s as defined in (j4.27j) below. 

Proposition 3.2. Let <j = q(e) = e a with a G (0, 1), and assume that 

• su P £ e(0,l) llXelloo < +oo. 

• w- lim £ _,o Xe^ 2n = u^ 2n G ^(M 2n ) 

• / R2 n Xe{w, q) {\w\ A + dist( ^g)2 ) dwdq<C < +oo. 

T/ien (|2.10p . (|2.1ip . and (|2.14p hold for the family of initial data {/0o,e} e e(o,i) • 

Proof. (|2.1ip follows from the first assumption and (|3.2p . 

Since ? = e a with a G (0, 1) we have that for all (w, q) G M 2n 

w- lim W £ ^ £ £> 2n = 6 {Wtq) , 

see |15[ Exemple III. 3], and so (|2.14p follows from our second assumption. 

To show that the third assumption implies (|2,10p . we notice that in this case (|2.10p can be 
written as follows 

/ Xe(w,q){H e ip £ w H £ \l) e w ) dwdq < +oo. (3.3) 
Since a < 1, and <f> G -£f 2 (K n ;C), by a simple computation we get 

(H £ ^ q ,H e ^ q ) < -(A x ^ q ,A x ^ q )+2(U^ q ,U^ q ) 

< C (l + \w\ 4 ) + C [ U(x) 2 -0 2 (?—^)dx. 

Since £/& is bounded, |[/ s ((/)| < C/dist(g, 5), and J* K „ |^(a;)| 2 dx = 1, a simple estimate analogous 
to the one in Section |4~41 shows that (|2.10p holds. We leave the details to the interested reader. □ 
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3.3 Conditions on the Wigner function 

Here we consider a general family of density matrices {po,e}ee(o,i) which satisfies the tightness 
conditions (|2 . 12[) and (|2. 13[) (so that f|2 . 14[) is satisfied up to the extraction of a subsequence). 
In the next proposition we show some simple sufficient conditions on the Wigner functions 
{We/9o,e}eG(o,i) m order to ensure the validity of assumptions (|2,10p and (|2.1ip . 

Proposition 3.3. Assume that 



max 



\ a \ t \/3\<w +1 \\d%d$W e p ,e\\oo <C< +00, 



• V (4- + U 2 (x) + \p\ 2 U(x) - zfAU(x)) W £ p ,e{x,p)dxd P < C < +00. 
Then (l2~T0l and (I2JTT) hold. 

Proof. Let us recall first that the Weyl symbol of an operator p of integral kernel p(x, y) is, by 
definition, given by 

o- £ (p)(x,p) := / p ( x+ y x -h e - i y^dy, 
that is equal to (27rs) n W £ p. Moreover, using (jA.3|) and (|A.4jl . it holds 

tr(p) = / W £ p(x,p)dxdp (3.4) 

Now, we remark that the first assumption gives (|2.1ip using Calderon-Vaillancourt Theorem 

Concerning (|2,10p . we will prove that 

sup tr(Hf/5o, £ ) < +00 
ee(0,i) 

(as observed in Section f2,31 this condition is slightly stronger than (j2.10|) ). To this aim, we first 
note that 

Hi = ^A 2 + U 2 --AU --UA. (3.5) 
6 4 2 2 v ' 

Moreover, let us observe that if p~\ and p~2 have kernels p\ and pi respectively, then the kernel 
associated to the operator p~\p~2 is given by J* pi(-,z)p2(z,-) dz. By this fact and (|3.4p . a simple 
computation shows that the identity 

ti(Ap e ) = / a E (A)(x,p)W E p , e (x,p)dxdp 

holds for any "suitable" operator A (here cr e (A) is the Weyl symbol of A). Hence, in our case, 

tr(#f Pe) = / o- e (H 2 )(x,p)W e p ,e(x,p)dxdp. 

JR 2n 
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We claim that the Weyl symbol of H 2 is 

a e (H 2 )(x,p) = + U\x) + \p\ 2 U(x) - ^AU(x). 

Indeed, let f(x,p) := \p\ 2 = a £ (—e 2 A)(x,p) and g(x,p) := U(x) = a e (U)(x,p). Then, using 
Moyal expansion, 

-4 ,-2 ,.2 

1 

= ^ X,P ^ +g(x,p) 2 + + 

where by definition 



<x £ (tf £ 2 )(x,p) = a £ (^ r A 2 + C/ 2 -^-AC/-^-C/A) (x,p) 



h^h 2 (x,p) := e i f^ 9 P '-^ 9 -')/ii(x,p)/i 2 (x / ,p / ) 



x'=x,p'- 



In our case, in the expansion of the exponential 



^( 9xdp ,-d pdx ,) = j2± (i £ -(d x d pl - d p d x , A 

j6N 1 



we can stop at the second order term, since f(x,p) = \p\ 2 . Therefore 



ng{x,p) = \p\ 2 U(x) - iep ■ VU(x) - ^-AU(x), 



2 



and 



g$f(x,p) = \p\ 2 U{x) + iep ■ VU(x) - V ^-AU(x). 



2 



This proves the claim and conclude the proof of the proposition. 



□ 



4 Proof of Theorem 12.1 



The proof of the theorem is split into several steps: first we show some basic estimates on the 
solutions, and we prove that the family W e p\ is tight in space and uniformly weakly continuous in 
time (this is the compactness part). Then we show that W e pf solves the Liouville equation (away 
from the singular set of the Coulomb potential) with an error term which converges to zero as e — > 
0. Combining this fact with some uniform decay estimate for W e p\ away from the singularity, we 
finally prove that any limit point is bounded and solves the Liouville equation. By the uniqueness 
of solution to the Liouville equation in the function space ([0, T]; L^IR 2 '™) n L°°(IR 2n )), we 
conclude the desired result. 

Let us observe that some of our estimates can be found [5] and [T]. However, the setting 
and the notation there are slightly different, and in some cases one would have to recheck the 
details of the proofs in [5j [T] to verify that everything works also in our case. Hence, for sake 
of completeness and in order to make this paper more accessible, we have decided to include all 
the details. 
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4.1 Basic estimates 

4.1.1 Conserved quantities 

The spectral decomposition of p\ is 

where 4>fl = e~ ttHe ^ £ (p^ solves (|1.2p . By standard results on the unitary propagator e ~ ttH rJ £ 
follows that 

E/4 E) <^4S> = E^ £ W £) > 

and 

E^W£¥=E/4 £ W?ii 2 (4-2) 

for all t £ R and e £ (0, 1). Therefore, using (|2,lL)p we have 

sup sup V ixf^H^f]) < +oc, (4.3) 
ee(o,i)*eRj^ 

sup sup E/4 e) ii^ff ii 2 < +°°- ( 4 - 4 ) 

£6(0,1) *eR^ 

4.1.2 A priori estimates 

From gH, (02]) and from the fact that U s > and C/ b £ L°°(R n ), follows that for all e € (0, 1) 
sup / C/ S 2 (x)pf(x,x)^< V/xf ||F £ #|| 2 + 2||C/ 6 || oc [y\ 

t£R JM. n f£ 

(4.5) 

supi J>f / \eV<j>f y {x)\ 2 dx <Y,$H<l>?H e <t>f) + ||C/ 6 ||oo. (4.6) 

Hence, by (j4 . 3[) and (|4.4p we obtain 

sup sup / U*{x)p £ t (x,x) dx < d (4.7) 
ee(o,i) teR JR™ 

and 

sup supV/iJ £) / leV^Cx)! 2 da; < C 2 . (4.8) 

ee(o,i)teR^ Jr™ 



(0f ,H e <t>f) + \\u b \v 



and 
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4.1.3 Propagation of ()2.11|) and consequences 

Observe that, by unitarity of e ttHs ^ e , we have, for all t £ R, 

—p e t < Cld. 



Hence, since 



w e pKv,p) = ^(^.ffi>. 



(see Appendix), using (j4.9|) we have 

Ce n C 
ee(o,i)teR ( 27r ) ( 27r ) 

(because ||$jj p || = £~ n / 2 ). Now, define for all x,y G M n and e, A > 

<fe, A)1 ,(s) = (v / 2e)" /2 (vrA)"/ 4 Gg(x - y). 

Observe that 

— (ge,\,y,Pt9eXy) = ^ \(9e,\,y, fifj \ 2 

therefore, since ||g e ,A,j/|| = x i by (|4.9p we have that 

2«/ 2 (vrAr/ 2 ^^Vg ) *Gg( y )| 2 <C. 



So 



sup sup sup XX IC< ) * G S(y)! 2 ^ » 
ee(o,i)te[o,r]yeR" ^ J ' A£ ^ ' 



4.2 Tightness in space 

Define = {y = (y\, . . . , y^) <E M. k : \yj\ < R, j = I, . . . , k}. We want to prove that 



lim sup sup / W e p £ (x,p) dxdp = 0. 

i?.^+oo E6(01) te ro.n jR2n\c ( „ 2n) 



ee(o,i) te[o,T] ./R 2 ™\C£ 
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Observe that for all R > 



sup sup 

ee(o,i) te[o,T] Jm. 2n \c^ n) 



W £ p\(x,p) dx dp < sup sup 



ee(0,l)*6[0,T] 



W £ pl(x,p) dx dp 



+ 



'x(R™\C<j n) ) 



W £ pl(x,p) dx dp 



so we can check the tightness property separately for the first and the second marginals of 
W £ p t e . From (|2,14p follows immediately that the family {W / e p £j o=§? 2n } £6 (o,i) is tight (because, by 
Prokhorov's Theorem, a family of nonnegative finite measures on R 2n is tight if and only if it is 
relatively compact in the duality with Ct,(M? n )). Therefore 



lim 



i \4 n) )xi 



W £ p £fl (x,p) dxdp = 0. 



(4.15) 



Let x e C(M n ), < x < 1 such that x(x) = if |x| < 1/2 and x{x) = 1 if \x\ > 1, and 
define Xr(x) ■= x{x/R). Observe that ||Vx/?||oo < C'/R and ||AxkIU < C'/R 2 . We define the 
following operator: 

4 £ V(z) = XR * G*?Xx)il>{x), $ G L 2 (R n ). 



Observe that 



and that [A§,H e ] = e 2 {A( XR * G? } )/2 + V( X R * G £ n) ) • V). So, using gS 



|tr(4 ) ^) = -Jtr([4 ) , J ff £ ]p*) 



d f 

= -77 XR(x)W £ p t £ (x,p)dxdp 

at J K 2n 



which gives 



/ 

J (i 



W £ pl(x,p) dx dp < 



< 



< 



XR(x)W £ pl{x,p) dxdp 

In 

XR(x)W £ po :£ (x,p) dxdp 

in 

W £ p 0i£ (x,p)dxdp + 

;n \C)j™ ) xK ™ 



R 2+ R 



R 2+ R 



T. 



Therefore, using (|4,15p . we get 



lim sup sup 

i?.^+oo ee (0,1) te[0,T] J(R n \C%g)xR n 



W £ pl(x,p) dxdp = 0, 



(4.16) 
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as desired. For the second marginal we observe first that 



\p\ 2 W £ pl(x,p) dx dp 



\p\ W £ p £ (x,p)dxdp+ — 



(4.17) 



and 



/ \p\ 2 W £ p t £ (x,p)dxdp = V/xf / 



E(e) 



therefore, using (|4.17p and (|4.8p . we have that 



\ fE 

(2vre) £ / 2 ^'* 



\eV(j) { -}{x)\ 2 dx 



n 



sup sup / \p\ W £ p £ (x,p) dx dp < C2 + 
ee(o,i)te[o,r] in 2 " 2 



and so 



1 



t e R k> := / ^(x.p)^ (x,p) 

JR 2 " 



is differentiable and 



sup sup 

ee(o,i) ten 



dt 



fedt) 



where Ca, is a constant depending only on (j). First observe that 



feA t )= W £ p £ t (x,p)<j> e (x,p)dxdp, 



where 



* G £ 2n \ Therefore, using (|2.2p . we have 
d 



dt 



feAt) 



<?e(Ub, P e t ){x,p)<t> £ {x,p) dx dp 
+ / ^s(U s , P £ t)(x,p)4> e (x,p) dx dp 

JR 2 " 

+ / (p-^x(t>e(x,p))W £ p £ t (x,p) dxdp. 

JR 2n 



dp 



(4.18) 



< sup sup / W £ pl(x,p) dx dp < ( C2 + ~ ) — >■ as i? — >■ +00. 

££(0,1) te[o,T] yR"x(R™\c^ n) ) ^ 

4.3 Weak Lipschitz continuity in time 

Here we prove that for all (j) € C^°(M? n ) the map 



(4.19) 



(4.20) 



(4.21) 



15 



For the first term it is easy to check that 
£e(U b , p £ t )(x,p)4> £ (x,p) dx dp 



< 



(27r) n 



\y\ sup \J r p (j) e \(x,y)dy. 



(4.22) 



.re li 



In the case of the Coulomb potential we follow a specific argument borrowed from [5j proof 
of Theorem l.l(ii)]), based on the inequality 



1 



1 



\z + w/2\ \z-w/2\ 



< 



w\ 



\z + w/2\\z - w/2\ 



(4.23) 



with z = (xi — xj) G R 3 , w = e{yi — yj) E R 3 . By estimating the difference quotients of U s as in 
(j4~23j) . using (|4~T|) we obtain 



^e(Us,pf)(x,p)(l)e(x,p) dxdp 



< C*d / \y\ sup ^(a/,!/)!^, (4.24) 



with C* depending only on the numbers Z\, . . . , an d Ci is the constant defined in (|4.7p . 
For the last term it is easy to see that 



(p ■ X7 x (f) £ (x,p))W £ p £ t (x,p) dxdp 



< 



1 



where 



(2vr) n j R n x . 
4>e{x,p) = p-V x (j) £ (x,p) 



sup \Tp4> £ \(x',y)dy, 



(4.25) 



Therefore we have only to bound 



\y\ sup \T p (j) E (x,y)\dy and / sup \T p <j> e \(x' ,y) dy 



with a constant depending only on 
For the first term 



|y| sup \F p 4> e (x,y)\ dy 



\y\ sup 



G^{x-^)J : p <l>{^,y)da! 



-V 2 e/4\ 



dy 



< / |y| sup |Jp</>(z,y)| < C^. 

il« zeK n 



For the second term 

sup \T p 4> e (x',y)\dy 



p y £ y*i> 5 y ) | t^y — i SUp 



dpdx'dp'e- ip - v (j){x',p')GW(p - p') (p ■ V x Gi n) {x - x') 

(4.26) 
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-ip-y 



Now observe that 

I dp dx'dp'e~ ip - y <l)(x',p')G^ (p - p') (p • V z G( n) (x - x') ' 
Jm. 3n \ J 

Tf dx'dp' d Xk G^(x-x')G^{p') f dp Pk cp{x\p-p') 

dx' (V x -T p g(x-x',y))G^(x') 
dx' {y-V x F p (p{x-x',y))Gf\x') 
where g(x,p) = p(j)(x,p). Now, since e £ (0, 1) 

( sup \T P Mx',y)\dy < [ sup ( dx' ( V • T p g(x - x', y))G^ {x') 

+ £ - [ sup [ dx' (y ■ VT^x - x', y))GW (x') 
< / dy sup |V ■ F P g(z,y)\ + t; / rf y 1 2/ 1 sup \V z T P 4>{z } y)\ 

JR n z£R. n * 7l» zeR n 



k=l 

-ey 2 /4 



+ - 



c 



(2) 



Therefore 



sup sup 

ee(o,i) *eR 



dt 



fed*) 



JJVC^JJoo (i) (1) 



(2) 



(2vr) n ' 



4.4 Uniform decay away from the singularity 

The singular set of U s is given by 

S= (J Sij, Sij = {x = (xi, ...,x M ,x)e (R 3 ) M x R n " 3M : Xi = xj for some i ^ j] , 

l<i<j<AI 

(4.27) 

and we have 



U a (x) > 



dist(x, S) ' 

where c > depending only on Z±, . . . , Zm- We want to prove that 

1 



lim sup 



N 4 + 



dist(x, S) 2 



W £ pl(x,p)dxdp < C. 



(4.28) 



(4.29) 
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We start with the second term: 



dx dp 



1 



dist(x, S) 2 



w £ p: 



Jb{ 



< 



1 

< - 

c 



») xR n dist(x,Dj z 
dx' 



dist(x',S) 2 
dx'U s (x') 2 pl(x',x') 



< 



where c is defined in (|4,28p . C\ is defined in (|4,7p . and we used (|4,28p . 
To prove the second estimate we observe that 

\p\ 4 W £ pl(x,p) dx dp < / \p\ A W £ p t e (x,p)dxdp 



+ 



\p\ 2 W £ pl. (x,p) dx dp + 



n(n + 2)e 2 



Thanks to (|4,18p . it suffices to control the first integral in the right hand side: 



/ \p\ 4 w £ p t £ (x,p)dxd P = y>f / 



dp 



(27re) £ / 2 ^'' i Ve 
£>f / \e 2 A^(x)\Ux 

H £( t>f t (x)\ 2 + U 2 {x)\^{x)\ 2 



dx, 



and the last term is uniformly bounded thanks to (|4,4p . ()4.7jl . and the boundedness of Ub- 

4.5 Limit continuity equation away from the singularities 

We want to prove that 



lim 

e-s>0 







/ 4>(x,p)W £ p £ (x,p) dx dp + (p(t) / b(x,p) ■ V(j)(x,p)W £ p £ (x,p) dx dp 



dt = 
(4.30) 

for all G C c °°(]R 2n \ (5 x IT)) and p G C c °°(0,T). Hence, recalling pU]) . we have to show that 

dxdp £' £ (U,p £ t )*G {2n) (x,p)<t>(x,p)+ [ dxdp VU(x)-V p (t>(x,p)W £ p t £ (x,p) = 0, 

Jr 2 ™ 

(4.31) 



lim sup 

e ^°te[o,T] 



and 



lim / dtip{t) [ dxdp^feV x -[W £ p £ t *Gi 2n) ](j)(x,p) = 0, (4.32) 
Jo JR 2 ™ 



for all G C, 



oo(m2n 



\ (S x M n )) and p G C£°(0,r). 
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4.5.1 Verification of (ELD) 

We can consider separately the contributions of Ub and U s . We start with the contribution of 
U s . We have to prove that 

lim sup / dxdp g £ (U s , p e t )*G ( l ? n) (x,p)4>(x,p) + dx dp VU s {x)-V p (f>(x,p)W £ p t £ {x,p) = 

(4.33) 



We know that 



for all (f) G C™(R 2n \ (S x R n )). 



lim sup 

£ ^°te[o,r] 



ip(x,p)W £ p t £ (x,p) dx dp — / tp(x , p)W £ pl(x , p) dx dp 



(4.34) 



for all y> G ^(IR 2 ™). 

First of all, we see that we can apply (|4.34p with (p(x,p) = X7U s (x) ■ V p 4>(x,p) to replace the 
integrals 

/ VU s (x) ■ V p <l)(x,p)W E tJj E dxdp 

with 

/ VU s (x) ■ V p (l)(x,p)W £ 'ifj £ dxdp 
in the verification of (|4.33p . Analogously, using (|4.7p and (|4.24p we see that we can replace 

f £e(U s , P i) * GjM(x,p)<j>(x, P ) dxdp 
jR 2n 

with 

/ & e (U s , pl)(x,p)(p(x,p) dxdp. 

jR 2n 

Thus, we are led to show the convergence 

lim sup / £ E (U s ,pl)<f) dxdp+ / VU s (x) ■ V p (p(x,p)W £ p t £ (x,p) dx dp = (4.35) 

for all G C~((M n \ S) x R n )). Since 



S e {U a ,p\)<f> dxdp 



U s (x + §?/) - U s (x - |y) t ( _ r 



R 2 ™ 



P £ ( x + y > x - y ) F P <i>(x, V) dxdy 



we can split the region of integration in two parts, where > 1 an d where \fe\y\ < 1. The 

contribution of the first region can be estimated as in (|4.24p . with 

C.I U^ M ,yU«( OH*)AM*, 

J{y/e\y\>l} x' JW- 
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which is infinitesimal, using (|4.7p again, as e — > 0. Since 

U a (x + %y) - U a (x - |y) 



VU a {x)-y 



uniformly as \fe\y\ < 1 and x belongs to a compact subset of W 1 \ S, the contribution of the 
second part is the same of 



(sy sy \ 

x + —,x - —) T p (j)(x, y) dxdy 



which coincides with 

/ W s (x) • V p 4>(x,p)W £ p t £ {x,p) dxdp. 
Now we consider the contribution of fT^. We have to prove that 



lim sup / £' £ (U b ,p t £ )(x,p)<t) £ {x,p) dxdp+ I VU b (x) ■ V p <p(x,p)W E p £ dx dp = (4.36) 
£ ^°te[o,T] . 



for all <f> E C c °°(M 2n ), where 
prove that 



* G £ 2n \ The proof of (|4,36p is divided in two parts: first we 



lim sup / £ £ (U h ,pl)(x,p)(f>{x,p) dxdp+ / VU b (x) ■ V p 4>(x , p)W £ p £ t dx dp = (4.37) 



for all 6 e C, 



oo(m2n\ 



and then, using the following estimate 



£e(Ub, p £ t )(x,p)(p(x,p) dx dp 



< 



(2vr)™ 



/ \y\ sup \T p (p\(x,y)dy. (4.38) 



for all if £ C^°(lR 2n ), we can replace (\> by (\> £ in the first summand of (|4.37p . obtaining (|4.36p . The 
proof of (|4.37p is achieved by a density argument. The first remark is that linear combinations 
of tensor functions (p(x,p) = 4>i(x)(j)2(p), with 0j 6 C£°(]R n ), are dense for the norm considered 
in (I4.38p . In this way, we are led to prove convergence in the case when 4>(x,p) = (j)i{x)4>2{p)- 
The second remark is that convergence surely holds if U b is of class C 2 (by the arguments in |15j . 
[5]). Hence, combining the two remarks and using the linearity of the error term with respect to 
the potential, we can prove convergence by a density argument, by approximating U b uniformly 
and in W 1,2 topology on the support of <j)\ by potentials V k 6 C 2 (M n ) with uniformly Lipschitz 
constants; then, setting A k = (U b — Vk)4>i and choosing a sequence in Lemma |4. II converging 
slowly to for k — > +oo, in such a way that 1 1 V^4^. 1 1 2 = o(A^ //4 ) for k — > +oo. In this way we 
obtain 

lim sup sup / S E (U h - V k ,p £ t )(x,p)(f)i(x)<l>2(p) dxdp = 0. 

fc->oo ee (o,l) te[0,T] JR 2 ™ 

As for the term in (|4.36p involving the Husimi transforms, we can use 04, lip to obtain that 

W £ p £ t V(U b (x) - V k {x)) ■ VMp)Mx) dxdp 

lim sup / \(pi(x)\\VU b (x) - VVjfc(x)| dx / |V</> 2 (p)| dp = 0. 

fe->oo JR™ JR™ 



lim sup sup sup 

fc->oo ee(0,l) te[o,T] 

C 



< 



(2tt) 
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So we need only to prove the following lemma: 

Lemma 4.1 (A priori estimate). For all A > ; we have that 



sup sup 

ee(o,i) te[o,T] 



£ e (U b , p e t )(x,p)4>i(x)4>2(p) dxdp 



(4.39) 



< 



>i i 



HV^Iloo sup \y\\My)- h *G { ^(y)\ + V\\\V AW^hh [ \u\G^\u) ^4.40) 



° llVAh f \z\\h\(z)dz + \\U b \U\VMoo [ \y\\h*G^\y)dy (4.41) 



(2vrA)™/ 4 

where A := i and C is i/ie constant in (|2.1ip . 

Proof. Set </>2 = Fp(f)2- Observe that since (|4.38p gives that 

sup sup [ £ e (U b ,p e t )Mx)Mp) dxdp- [ £ £ (U b ,p E t )Mx)Mp)e~ lpl2X dxdp 
ee(o,i)te[o,T] Jm 2 ™ Jm 2 " 

< ll^illillV^Hoo sup \y\\My) - h * G[ n) (y)\ 



we recognize the first error term in (j4J39]). So we have only to estimate 

£ £ (U b ,p £ t )Mx)Mp)e- lpl2X dxdp 



sup sup 

£6(0,1) t6[0,T] 



(4.42) 



Observe that 

I g £ (U b ,p e t )Mx)Mp)e~ lpl2X dxdp = I e>t + II e>t - III £>t 

JR 2n 

where 

hy-=\ + IV) - MX - %V) 4>, * Gf (y)pt(x + £ -y,x- £ -y) dxdyd, (4.43) 

He,f.= [ U b (x + £ -y) Mx) ~ MX + § ^ 2 * Gt\y) P l (x + £ -y,x- £ -y) dxdy, (4.44) 

JM. 2n Z £ Z Z 

UI £ y.= -\ U b (x - £ -y) Mx) ~ MX ~ ^ y) h * G^(y)pt(x + £ -y,x- £ -y)dxdy. (4.45) 

jR2n Z £ Z Z 



Observe first that 



sup sup \II £>t \ + \III £jt \ < ||C/ 6 |UI|V0i|U / \y\\4>2*G { ^\(y)dy. 

£6(0,1) t6[0,T] JR" 

The estimate of I £t t is more delicate: we first perform some manipulations of this expression, then 
we estimate the resulting terms with the help of (14. 13ft . 
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We expand the convolution product and make the change of variables 



e e 

u = x H — y v = x y 

2 2 



to get 



dudvdz 



A(u) - A(v) \ez-(u-v)\ 



A p 6 t (u,v)(j)2{z) 



£ ' 1 (7tA)"/ 2 £ 

= lH^f I (A$)*G^v + ez)l^)4> 2 (z)dvdz 



IE 



(e) 



* G?g)(u + ez)<i>ll{v)<h{z) dvdz 



J- v-^ 



(e) 



(^) * G<£(u + ez) - A(v + ez)(45 * G^)(v + ez) <^}{v)fc{z) dvdz 



/xf / [A(t; + ez) - (0$ * G^' 2 )(v + ez)0fj(v)Mz) dvdz. 



(e) _ ^(n) 



(4.46) 



Now let us estimate the first summand in (|4.46|) 



1 

j'eN 



(6) 



(^g ) * eg (« + ez) - A(v + ez) (</>f> * ) (u + ez) <jlf> (v) fa (z) cfccfe 



(e) ^(n) • 



<iz0 2 (z) / dudt, + ^ - ") - ^ + £ ^ G 0O w ff W g {V + £Z _ u) 



< IIVAIU / dz\fa{z)\ I dudv^GW(u)\<l>fl(v)\\(l)f>(v + ez-u)\ 

< x/AllVAIUH^IIi / lulG^H^. 

For the second summand in (j4.46j) . using (I4.13p . we have 

- J2 »f [ [A(v + ez) - * )(« + ez)$}{v)fc{z) dvdz 



< 



< 



E". 



R 2n 

A(w + ez) - 



| (^5 * Gg )(u + «0|| 45 («) 1 1 02 (*) | dvdz 



(«), 



C 



||VA|| 2 / \zU 2 (z)\dz. 



This completes the estimate of the term in (|4,43p and the proof. 



□ 
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4.5.2 Verification of (jiT32j) 

This is easy, taking into account the fact that 

W £ pl * Gf n \x,p) ■ V x <j)(x,p) dxdp= [ W £ plV x ■ [4> * Gf n \ dxdp 

are uniformly bounded (recall that G e 2n \ defined in (|2.7p . are probability densities). 
4.6 Proof of Theorem I2TT1 

Define : [0,T] -> ^(R 2n ) as W t [e) := W s p\££ 2n for all e G (0,1) and t G [0,T]. Using 

(|4.14p . (|4,19p and Ascoli-Arzela Theorem, one can prove easily that there exist a subsequence 
{ W{£k) ) km and W ■ [°' T ] -> ^(K 2n ) such that 

lim sup d^(W t {£k) ,W t ) = 0. (4.47) 

We now prove the following assertions: 

(i) W : [0,T] ^(]R 2n ) is weakly continuous and, for all t G [0,T], W t = #^ 2ri for 
some function W t G L^M 2 ™) n L°°(R 2n ). Moreover W t > and sup 4e[0iT] ||^||ii( K 2«) + 
||#t||L-(R2n ) <C. In particular, W G £+([0, T]; L 1 (R 2n ) n L°°(]R 2n )). 

(ii) b G Lj^^O,! 1 ) x ]R 2n ; , so the continuity equation (|2.8p with ut = Wt makes sense; 

(iii) W solves (|2.8p in the sense of distributions on [0,T] x R 2n ; 

(iv) For any <j> G C£°(M 2n ), i ^ / R2 „ <£dWt belongs to ^([O.T]). 

Proof of (i): Observe that (|4.47p implies that : [0, T] — )• ^(lR 2n ) is weakly continuous because 
it is uniform limit of the weakly continuous maps W^ £k \ The second part of the proposition 
follows immediately from (|4,lip . Indeed, for all <fi G L 1 (BL 2n ), 

sup sup f (j)(x,p)W £ pf(x,p) dxdp< f 4>(x,p) dxdp (4.48) 
e6(o,i) te[o,T] Jm 2n (2iT) n J R 2 n 

and so 

sup / <f>(x,p)dWt(x,p) < (f)(x,p) dxdp. (4.49) 
te[o,T] Jm 2 ™ ( 27r ) Jr 2 " 

Proof of (ii): The estimate b G I^ c ((0, T) x M 2n ; dWidt) follows easily from (ITOD and (ITOD . 

Proof of (iii): First we prove that W solves (|2.8p in R 2n \(Sx M n ), where S is the singular 
set of U s defined in (|4.27p . Unfortunately this does not follow immediately by (|4.30p because we 
have no information about the singular set £ of VC/ft, so we cannot control the limit k — > oo of 

T r 

dt<p(t) / dxd P VU b (x) ■ V p <i>{x,p)W eh p\ {x,p) 



o 



23 



in (|2,8p with (|4,47p . But we can proceed by a density argument because, using the regularity 
conditions (j4.48[) and (|4,49p . we can approximate VJ7f, in L 1 on suppcfi by bounded continuous 
functions. 

In order to prove that # solves (JUJ) in [0, T] x R 2n we use (1429]) to obtain that 

SU P / i ■ TT dW t (x,p)dt < +oo. (4.50) 
te[o,T] d\st(x, S) 2 

Observe that (1PU1) implies that W t {S x R n ) = for every f £ (0,T). The proof of the global 
validity of the continuity equation uses the classical argument of removing the singularity by 
multiplying any test function (f> G C£°(K 2n ) by Xki where Xk( x ) = x(fcdist(a;, S)) and x is a 
smooth cut-off function equal to on [0, 1] and equal to 1 on [2, +oo), with < x' < 2. If we 
use </>Xfc as a test function, since Xk depends on x only, we can use the particular structure of b, 
namely b(x,p) = (p, — VC/(x)), to write the term depending on the derivatives of Xk as 

k / <px'{kdist{x,S))p ■ Vdist(x, S)dW t (x,p)dt. 

If K is the support of 0, the integral above can be bounded by 

2 sup / kdW t (x,p)dt< 8maX « lp(pl [ 1 c , dW t (x,p), 

K J{x£K:kdist{x,S)<2} k JK dlSt {X,b) 

and the right hand side is infinitesimal (uniformly in t) as k — > oo. 

Proof of (iv): Since the distributional derivative of £ i — >- J^ 2n (ftWfdxdp is given by J" R2 n b ■ 
VcpdWt, we have to show that the map 

[ b-V(j)dW t 

jR 2n 

is continuous. Observing that the map t t-> Wt is weakly continuous and Wt = Wt^£ 2n with 
# G ([0, T];L 1 (R 2n ) n L°°(IR 2n )), the only delicate term is 

( VU s (x)-V p $(x,p)dW t . 

Define the nonnegative Hamiltonian function T~L = \p\ 2 /2 + U +\\Ub\\oo- Taking the limit in ([4.29)1 
as e —> we easily deduce that 

sup f U 2 dW t < C sup I (l + |p| 4 + U 2 (x)) dW t < +oo. 

Since the Hamiltonian is preserved by the Liouville dynamics (under our assumptions on the 
potential, this fact is contained in the proof of [TJ Theorem 6.1]), the above bound implies 

sup / H 2 dW t = [ n 2 dW as N -»• oo. 

te[0,T] J{H>N} J{H>N} 
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As U s < H, this implies 

sup f U 2 dW t < I H 2 dW ^ as N -> oo. 

te[0,T] J{Us>N} J{H>N} 

Hence, if we define the sets An '■= {U s < N}, the functions 

t^f N (t) := [ VU s (x)-V p ct>dW t 
J A N 

are continuous and converge uniformly to JJ^n W s {x) ■ W p (pdWt as N — > oo. This proves (iv). 

To conclude the proof of the theorem, recalling that W denote the unique distributional 
solution of JZaj in L+ ([0, T];L 1 (R 2n ) n L°°(E 2n )) starting from oj^ 2n (see (TJ Theorem 6.1]), 
we have proved W = W , and so 

lim sup d&(W ek pl, & 2n ,W t ^ 2n ) = 0. (4.51) 

Since the limit #jJz? 2n is independent of the chosen subsequence, this implies the convergence of 
the whole family, namely 

lim sup d&{W £ p l £ ^ 2n ,W t ^ 2n ) = 0, (4.52) 
£ ^°te[o,T] 

as desired. 

A Notations and some notions about density operators 

A density operator on L 2 (R n ) is a positive, self-adjoint, trace-class operator, namely p = p* , p > 
and tr(p) = 1, where the trace is defined as follows: 

: =X^''W) ( A -!) 

with {^}j g N is an y orthonormal basis of L 2 (R n ). It can be shown that each density operator p 
is a compact operator, so it can be decomposed as follows 

p = J2 X Mi-)iPj (A.2) 

where < Ai < A2 < • • • < 1, and {V'jljeN is a orthonormal basis of eigenvectors of p. Therefore 
p is an integral operator and its kernel is 
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so that 

pij)(x) = / p(x,y)ip{y) dy. 
Observe that the trace condition on p can be expressed as follows in terms of its kernel 

tr(p) = / p(x,x)dx = l. (A.3) 



The Wigner transform of p is defined as 

W £ p(x,p) := - — — / p(x + -y, x — -y)e~ tpy dy, 
and the Husimi transform of p as 



W £ p := W £ p * Gf n \ G? n \x,p) := Gt\x)G { p\p) 



e 



(vre) n 

It is easy to check that the marginals of W e p are 

/ W e p(x,p) dp = p(x,x) and / W e p(x,p) dx = . F (-,-) (A.4) 

where 

J r p(q,q)= [ p(u,u)e- iq - u du. (A.5) 



(A.6) 



Similarly the marginals of W £ p are 

f W £ p{x,p)dp= I' p(x-x',x-x')G ( i l \x')dx' 



and 



/ W £ p(x,p) dx = -i— / Tp P —f) Gi n \p') dp'. (A.7) 

Moreover, the Husimi transform is nonnegative: indeed (see for instance |15j). 

W^(x,p) = ^\(i;,f X)P )\ 2 , (A.8) 

where (•, •) is the scalar product on L 2 (M ra ) and 

■= jjy^ e ~ lx ~ yl2/(2£)e ~ l(p ' y)/e G L2 ( Rn )' WtiJi = L 

Hence W e ip > 0, and using the spectral decomposition (jA.2|) one obtains the non-negativity of 
W £ p for any trace-class operator p. Moreover, combining fjA.3|) and (1A.6|) . it follows that W £ p is 
a probability measure. 
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